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Abstract 

Mathematical  modelling  is  an  essential  tool  for  the  design  of  solid  oxide  fuel  cells  (SOFCs).  The  present  paper  aims  to  report  on  the 
development  of  a  dynamic  anode-supported  intermediate  temperature  direct  internal  reforming  planar  solid  oxide  fuel  cell  stack  model,  that 
allows  for  both  co-flow  and  counter-flow  operation.  The  developed  model  consists  of  mass  and  energy  balances,  and  an  electrochemical 
model  that  relates  the  fuel  and  air  gas  composition  and  temperature  to  voltage,  current  density,  and  other  relevant  fuel  cell  variables.  The 
electrochemical  performance  of  the  cell  is  analysed  for  several  temperatures  and  fuel  utilisations,  by  means  of  the  voltage  and  power  density 
versus  current  density  curves.  The  steady-state  performance  of  the  cell  and  the  impact  of  changes  in  fuel  and  air  inlet  temperatures,  fuel 
utilisation,  average  current  density,  and  flow  configuration  are  studied.  For  a  co-flow  SOFC  operating  on  a  10%  pre-reformed  methane  fuel 
mixture  with  75%  fuel  utilisation,  inlet  fuel  and  air  temperatures  of  1023  K,  average  current  density  of  0.5  A  cm-2,  and  an  air  ratio  of  8.5,  an 
output  voltage  of  0.66  V  with  a  power  density  of  0.33  W  cm'2  and  a  fuel  efficiency  of  47%,  are  predicted.  It  was  found  that  cathode  activation 
overpotentials  represent  the  major  source  of  voltage  loss,  followed  by  anode  activation  overpotentials  and  ohmic  losses.  For  the  same  operating 
conditions,  SOFC  operation  under  counter-flow  of  the  fuel  and  air  gas  streams  has  been  shown  to  lead  to  steep  temperature  gradients  and 
uneven  current  density  distributions. 
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1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  are  energy  conversion 
devices  that  produce  electricity  and  heat  directly  from  a 
gaseous  or  gasified  fuel  by  electrochemical  combination  of 
that  fuel  with  an  oxidant.  A  SOFC  consists  of  an  intercon¬ 
nect  structure  and  a  three-layer  region  composed  of  two  ce¬ 
ramic  electrodes,  anode  and  cathode,  separated  by  a  dense 
ceramic  electrolyte  (often  referred  to  as  the  PEN — Positive- 
electrode/Electrolyte/Negative-electrode).  SOFCs  operate  at 
high  temperatures  and  atmospheric  or  elevated  pressures, 
and  can  use  hydrogen,  carbon  monoxide,  and  hydrocar¬ 
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bons  as  fuel,  and  air  (or  oxygen)  as  oxidant.  In  the  cell, 
the  oxygen  ions  formed  at  the  cathode  migrate  through  the 
ion-conducting  electrolyte  to  the  anode/electrolyte  interface 
where  they  react  with  the  hydrogen  and  carbon  monoxide 
contained  in  (and/or  produced  by)  the  fuel,  producing  wa¬ 
ter  and  carbon  dioxide  while  releasing  electrons  that  flow 
via  an  external  circuit  to  the  cathode/electrolyte  interface  [1]. 
An  oxygen  chemical  potential  gradient  is  formed  across  the 
electrolyte  by  burning  fuel  with  oxygen  on  the  anode,  while 
leaving  the  cathode  exposed  to  abundant  atmospheric  oxy¬ 
gen. 

While  conventional  high-temperature  SOFCs  generally 
operate  between  1073  and  1273  K,  a  number  of  research 
groups  are  presently  focusing  on  intermediate  temperature- 
solid  oxide  fuel  cells  (IT-SOFCs).  These  typically  operate 


P.  Aguiar  et  al.  /  Journal  of  Power  Sources  138  (2004)  120-136 


121 


Nomenclature 

ASR 

area  specific  resistance  (£2  m2) 

cp 

specihc  heat  capacity  of  the  gas  streams,  PEN 

structure,  and  interconnect  (kJkg-1  K-1) 

Ci 

molar  concentration  of  component  i  (mol  m-3) 

dfo 

hydraulic  diameter  of  a  rectangular  cross- 

section  gas  channel  (m) 

^bulk,  i 

bulk  diffusivity  coefficient  (m2  s-1) 

Deff, electrode  electrode  effective  diffusivity  coefficient 

(m2  s-  ) 

Z^Knudsen,  i  Knudsen  diffusivity  coefficient  (m2  s  1 ) 

Anolecuiar,  i  molecular  diffusivity  coefficient  for  com- 

ponent  i  (m2  s_1) 

Dx 

axial  diffusion  coefficient  (m2  s-1) 

£a 

activation  energy  (kJ  mol-1) 

^electrode 

activation  energy  of  the  exchange  current  den- 

sity  (kJ  mol-1) 

F 

Faraday’s  constant  (C  mol-1) 

Fair 

molar  flow  rate  of  the  air  stream  (mol  s-1) 

^fuel 

molar  flow  rate  of  the  fuel  stream  (mol  s-1) 

K 

air  channel  height  (m) 

fuel  channel  height  (m) 

j,  j 

average  and  local  current  density  (Am-2) 

yo, electrode  exchange  current  density  (A  m  ~) 

k-d,b  ka,PEN  air  channel  heat  transfer  coefficient 

(kJm-2  s-1  K-1)  (Table  4) 

^electrode 

pre-exponential  factor  of  the  exchange  current 

density  (A  m-2) 

kfj,  k\\ 

pen  fuel  channel  heat  transfer  coefficient 

(kJm-2  s-1  K-1)  (Table  4) 

ka 

equilibrium  constant  for  hydrogen  adsorption 

(Table  3) 

ks 

equilibrium  constant  for  steam  adsorption 

(Table  3) 

ko 

pre-exponential  constant  (Table  3) 

L 

system  length  (m) 

LHV 

lower  heating  value 

n 

number  of  electrons  participating  in  the  elec- 

trochemical  reaction 

Nu 

Nusselt  number 

P 

pressure  (bar) 

Pe 

Peclet  number  Pe  =  uL/Dx 

Pi 

partial  pressure  of  component  i  in  relevant  gas 

channel  (bar) 

Pi.  TPB 

partial  pressure  of  component i  at  relevant  three 

phase  boundary  (bar) 

Rsofc 

power  density  (W  m-2) 

Rk 

rate  of  reaction  k  (mol  m-2  s-1) 

Rohm 

total  cell  resistance,  including  both  ionic  and 

electronic  resistances  (f2m2) 

m 

gas  constant  (kJ  mol-1  K-1) 

t 

time  (s) 

T 

temperature  (K) 

Ma 

air  velocity  (ms-1) 

Uf 

fuel  velocity  (ms-1) 

u 

potential  (V) 

u°c  p 
< 

open-circuit  potential  (OCP)  (V) 

OCP  for  standard  temperature  and  pressure  and 
pure  reactants  for  the  Pb  oxidation  reaction  (V) 

^fuel 

fuel  utilisation  factor 

w 

system  width  (m) 

X 

axial  coordinate  (m) 

yi 

molar  fraction  of  component  i 

y? 

molar  fraction  of  component  i  at  inlet 

Greek  letters 

a 

transfer  coefficient 

fiii  Pi 

methane  and  steam  order  of  reaction  (Table  3) 

Pi 

number  of  surface  sites  required  for  methane 
adsorption  (Table  3) 

(—AH)k 

enthalpy  change  of  reaction  k  (kJ  mol- 1 ) 

e 

emissivity 

k-dir 

air  ratio 

kd,  kf , 

kpEN-  k\  thermal  conductivity  of  the  gas 
streams,  PEN  structure,  and  interconnect 
(kJ  m-1  s-1  K-1) 

Sp 

electrode  porosity 

esoFC 

fuel  cell  efficiency 

^?act 

activation  overpotential  losses  (V) 

Oconc 

total  concentration  overpotential  losses  (V) 

^conc, electrode  concentration  overpotential  losses  (V) 

a 

Stefan-Boltzmann  constant  (W  m  2  K  4) 

a anode  -  ^cathode  electronic  conductivity  of  the  anode  and 

cathode  (£2  m 

^electrolyte  ionic  conductivity  of  the  electrolyte 

(f2-1  m-1) 

Pa,  Pi,  Ppen,  Pi  density  of  the  gas  streams,  PEN  struc- 

ture,  and  interconnect  (kg  m  3) 

^anode,  rcathode  thickness  (m) 

^electrolyte »  "^PEN?  ^1?  ^tortuosity  tortuosity 

Vik 

stoichiometric  coefficient  of  component ;  in  re¬ 
action  k 

Superscripts 

0 

feed  conditions  (fuel  and  air  channel  inlet) 

Subscripts 

a 

air  channel 

act 

activation 

cone 

concentration 

f 

i 

fuel  channel 
component 

i 

interconnect 

k 

reaction 

Ohm 

Ohmic 

PEN 

PEN  structure 

TPB 

three-phase  boundary 
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between  823  and  1073  K,  allowing  for  a  wider  range  of 
materials  and  more  cost-effective  SOFC  fabrication,  par¬ 
ticularly  in  relation  to  the  interconnections  and  balance  of 
plant  (BoP).  High-temperature  SOFCs  are  generally  all  ce¬ 
ramic  while  intermediate  temperature  SOFCs  are  metal- 
ceramic  and  use  stainless  steel  interconnects  instead  of  more 
expensive  high-chrome  alloys  or  oxides.  A  typical  SOFC 
electrolyte  is  yttria-stabilised  zirconia  (YSZ),  an  oxide  ion 
conductor  at  elevated  temperatures.  The  anode  is  usually  a 
nickel/zirconia  cermet,  which  provides  high  electrochem¬ 
ical  performance,  good  chemical  stability,  and  low  cost, 
and  the  cathode  a  perovskite  material,  such  as  strontium 
doped  lanthanum  manganite,  often  mixed  with  YSZ  in  the 
form  of  a  composite  [2,3].  Two  main  SOFC  designs  are  the 
electrolyte-supported  and  electrode-supported  SOFC.  In  the 
former  design,  the  electrolyte  is  the  thickest  component  and 
works  as  the  support  structure.  Electrolyte-supported  cells  are 
only  suitable  for  high-temperature  operation  where  the,  of¬ 
ten  large,  ohmic  losses  can  be  reduced.  Electrode-supported 
SOFCs  have  been  developed  in  an  attempt  to  minimise  such 
ohmic  losses  under  lower  temperature  operation,  as  in  the 
case  of  IT-SOFCs.  In  these  cells,  one  of  the  two  electrodes 
is  the  thickest  component  and  support  structure,  while  the 
electrolyte  is  required  to  have  high  ionic  conductivity  and/or 
small  thickness.  However,  it  is  usually  observed  that,  despite 
a  low  ohmic  contribution,  the  area  specific  resistance  of  IT- 
SOFCs  may  be  larger  than  high-temperature  SOFCs.  This 
is  due  to  activation,  and  perhaps  concentration  overpoten¬ 
tials,  that  can  often  outweigh  the  ohmic  contribution.  Thus, 
many  of  the  recent  efforts  in  fuel  cell  technology  develop¬ 
ment  have  been  devoted  to  reducing  the  thickness  of  critical 
cell  components  while  refining  and  improving  the  electrode 
structure  and  electrolyte  phase,  with  the  aim  of  obtaining  a 
higher  and  more  stable  electrochemical  performance,  while 
lowering  cost  [4-6]. 

Most  fuel  cells  need  to  convert  a  hydrocarbon  primary 
fuel  into  a  hydrogen-rich  gas  required  for  the  electrochem¬ 
ical  reaction  on  the  anode  side.  One  option  for  converting 
the  hydrocarbon  feed  is  to  do  it  indirectly  in  a  fuel  process¬ 
ing  system,  such  as  an  external  catalytic  steam  reformer  or 
partial  oxidation  reactor  [7,8].  In  the  most  common  case  of 
a  reformer,  heat  needs  to  be  available  to  drive  the  steam  re¬ 
forming  reaction.  One  method  of  achieving  this  is  to  feed 
the  exhaust  gases  from  the  anode  and  cathode  into  a  burner 
where  the  excess  fuel  is  combusted.  The  heat  generated  in 
the  burner  can  then  be  used  to  preheat  both  the  steam  and  the 
fuel,  and  provide  the  heat  needed  in  the  reformer.  However, 
a  more  elegant  and  efficient  method  of  providing  heat  for  re¬ 
forming  in  fuel  cell  systems  is  to  carry  out  the  reforming  in 
the  cell  stack  [9,10].  Catalytic  steam  reforming  of  hydrocar¬ 
bons  is  typically  performed  between  1023  and  1 173  K,  and  is 
thus  compatible  with  SOFCs.  These  operate  at  high  enough 
temperatures  for  the  endothermic  steam  reforming  reaction 
to  be  carried  out  within  the  stack,  giving  rise  to  an  internal 
reforming  SOFC  (IR-SOFC).  Depending  on  the  operating 
conditions,  the  heat  for  the  reforming  reaction  can  vary  from 


40  to  70%  of  the  total  heat  produced  in  a  fuel  cell.  Thus, 
internal  reforming  fuel  cells  eliminate  the  requirement  for  a 
separate  fuel  reformer,  leading  to  a  more  attractive  and  effi¬ 
cient  system  design;  the  requirement  for  cell  cooling,  which 
is  usually  achieved  by  flowing  excess  air  through  the  cath¬ 
ode,  is  also  significantly  reduced  [11,12].  There  are  two  main 
approaches  to  internal  reforming  within  a  fuel  cell:  indirect 
or  integrated  (HR)  and  direct  (DIR)  internal  reforming.  In  the 
first  approach,  the  reformer  section  is  separate  but  adjacent  to 
the  fuel  cell  anode  and  in  close  thermal  contact  with  it.  In  the 
latter  approach,  methane  is  fed  directly  into  the  cell  and  the 
reforming  takes  place  directly  on  the  anode  [9].  One  advan¬ 
tage  of  HR  is  that  the  reformer  and  cell  environments  do  not 
have  a  direct  physical  effect  on  each  other;  however,  the  con¬ 
version  of  methane  to  hydrogen  is  not  promoted  to  the  same 
extent  as  with  direct  internal  reforming.  For  a  DIR  configu¬ 
ration,  part  of  the  steam  required  for  the  reforming  reaction 
can  be  obtained  from  the  fuel  cell  electrochemical  oxidation 
of  hydrogen  and,  because  of  the  continuing  consumption  of 
hydrogen,  the  equilibrium  of  the  reforming  reaction  may  be 
further  shifted  to  the  right,  increasing  the  methane  conversion 

[13]  and  leading  to  a  more  evenly  distributed  load  of  hydro¬ 
gen.  However,  the  DIR  approach  requires  an  anode  material 
that  possesses  good  catalytic  properties  for  steam  reforming 
as  well  as  for  effective  generation  of  power  [12].  In  a  SOFC, 
the  reforming  may  occur  internally:  it  has  been  shown  that 
the  state-of-the-art  SOFC  nickel/zirconia  cermet  anodes  can 
provide  sufficient  activity  for  the  steam  reforming  and  shift 
reactions  without  the  need  for  any  additional  catalyst  [7,9]. 
A  problem  related  to  DIR  is  carbon  deposition  on  the  anode 
and  subsequent  electrocatalyst  deactivation,  leading  to  loss  of 
cell  performance  and  poor  durability.  The  high  steam/carbon 
ratio  typically  used  in  conventional  steam  reformers  to  sup¬ 
press  carbon  formation  is  unattractive  as  it  lowers  the  elec¬ 
trical  efficiency  of  the  fuel  cell  by  steam  dilution  of  the  fuel 

[14] .  Hence,  advanced  anode  materials  that  allow  for  direct 
internal  reforming  at  low  steam/carbon  ratios  could  offer  sig¬ 
nificant  benefits.  A  second  problem,  mostly  related  to  high- 
temperature  SOFCs,  is  the  strong  cooling  effect  caused  by 
the  highly  endothermic  reforming  reaction  that  can  gener¬ 
ate  large  temperature  gradients  across  the  cell  and  limit  the 
amount  of  internal  reforming  allowed  in  practice.  It  has  been 
shown  that  in  high-temperature  IR-SOFCs,  all  the  methane 
is  usually  completely  reformed  within  a  small  distance  from 
the  anode  entrance  [15-18].  IT-SOFCs  also  offer  significant 
advantages  in  terms  of  internally  reforming  methane,  as  the 
lower  temperature  naturally  reduces  the  reforming  reaction 
rate,  without  the  need  for  any  catalyst  modification  or  com¬ 
plex  cell  design. 

Mathematical  modelling  is  an  essential  tool  in  the  design 
of  fuel  cell  systems,  as  it  is  important  to  understand  the  re¬ 
sponse  of  a  cell  stack  under  normal  or  transient  conditions. 
Such  models  can  provide  a  picture  of  cell  stress,  potential, 
current  density,  and  temperature  as  functions  of  position  and 
time,  for  various  cell  configurations  and  operating  conditions, 
and  be  used  to  examine  the  effects  of  changes  on  one  or 
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more  variables  and  the  relative  system  sensitivity  to  relevant 
design  parameters.  A  model  is  also  useful  in  predicting  the 
effects  of  altering  process  variables  and  in  using  that  infor¬ 
mation  to  optimise  cell  performance.  The  main  aims  of  this 
paper  are  to:  present  a  dynamic  model  of  a  planar  anode- 
supported  intermediate  temperature  SOFC  stack  with  direct 
internal  reforming;  investigate  the  steady-state  performance 
of  such  SOFC  system;  and  analyse  the  impact  of  changes 
in  inlet  temperature,  fuel  utilisation,  average  current  density, 
and  flow  configuration.  Although,  the  full  dynamic  SOFC 
model  has  been  developed  and  is  explained  here,  the  present 
paper  only  explores  the  cell  steady-state  performance.  A  sub¬ 
sequent  publication  will  deal  with  its  dynamic  aspects  [19]. 
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Fig.  1 .  Schematic  side  view  of  a  co-flow  planar  SOFC  stack. 


2.  SOFC  mathematical  model 

There  have  been  several  publications  focusing  on  mod¬ 
elling  the  performance  of  high-temperature  SOFCs.  Such 
models  can  be  for  different  geometries  (tubular,  planar, 
or  monolithic),  PEN  structures  (electrolyte  or  electrode- 
supported),  and  flow  configurations  (cross-,  co-,  or  counter¬ 
flow);  can  range  from  one  to  three-dimensional;  and  can 
consider  or  ignore  internal  reforming  or  various  other  phe¬ 
nomena  [1,15,20-25].  Even  though  most  of  the  work  found 
is  for  SOFC  steady-state  performance,  dynamic  issues  have 
also  been  addressed  by  some  authors  [26-29].  The  dynamic 
intermediate  temperature  SOFC  model  developed  here  is  one¬ 
dimensional  and  allows  for  co-  and  counter-flow  operation, 
although  the  equations  presented  next  are  for  the  co-flow  case 
only. 

To  produce  a  useful  voltage,  a  complete  SOFC  system  can 
not  be  made  up  of  a  single  solid  oxide  fuel  cell,  as  its  voltage 
is  quite  small  (about  0.7  V  when  drawing  a  useful  current), 
but  must  consist  of  several  repeating  electrochemical  cells  in 
a  module,  connected  both  in  series  and/or  in  parallel  and  as¬ 
sembled  to  compose  a  stack.  However,  models  of  SOFCs  are 
usually  developed  for  the  smallest  unit  cell  or  module.  Such 
models  can  then  be  used  to  compute  the  operating  conditions 
of  the  whole  stack  provided  that  the  boundary  conditions  for 
each  unit  cell  are  independent  and  that  adequate  boundary 
conditions  are  established  for  the  stack  exterior.  In  the  case 
reported  here,  the  repeating  single-cell  is  considered  to  be 
in  the  centre  of  a  large  stack,  such  that  no  end  effects  are 
present,  and  to  be  formed  when  two  interconnect  plates  are 
placed  above  and  below  an  electrochemically  active  tri-layer 
cell  composed  of  the  anode,  electrolyte,  and  cathode  (PEN 
structure).  Interconnect  plates  normally  provide  the  internal 
gas  flow  manifolds  and  flow  passages  above  and  below  the 
PEN;  the  individual  gas  flow  channels  are  neglected  in  the 
present  model.  Fig.  1  shows  a  schematic  view  of  a  vertical 
slice  through  a  y-z  plane  of  a  multiple-cell,  planar,  co-flow 
stack,  where  the  fuel  and  air  gas  streams  flow  along  the  x- 
direction.  For  the  purpose  of  modelling,  it  is  considered  that 
the  SOFC  is  composed  of  the  fuel  and  air  channels,  the  PEN 
structure,  and  the  interconnect.  The  model  developed  con¬ 


sists  of:  two  mass  balances,  to  the  fuel  and  air  channels;  four 
energy  balances,  to  the  fuel  and  air  channels,  PEN,  and  inter¬ 
connect;  and  an  electrochemical  model  that  relates  the  fuel 
and  air  gas  composition  and  temperature  to  the  voltage,  cur¬ 
rent  density,  and  other  cell  variables.  Pressure  drop  along  the 
channels  is  neglected.  Note  that  for  the  development  of  such 
a  model,  several  simplifying  assumptions  (namely,  the  vari¬ 
ation  of  some  gas/solid  properties  along  the  fuel  cell  length 
and  with  time)  were  made  to  obtain  a  more  computationally 
tractable  dynamic  model.  Although  it  is  believed  that  any  dif¬ 
ferences  between  the  results  predicted  with  this  model  or  with 
a  more  refined  one  will  not  be  significant,  such  a  model  will 
also  be  developed  in  the  near  future  and  the  results  obtained 
reported. 

2.1.  Mass  balances 

For  the  mass  balance  in  the  fuel  channel,  the  chemical 
species  considered  are  CH4,  HiO,  CO,  H2,  and  COt,  while 
for  the  air  channel  the  chemical  species  are  O2  and  N2  (as 
the  oxidant  is  considered  to  be  air  and  not  pure  oxygen).  The 
molar  flux  in  these  gas  channels  is  considered  convective  in 
the  flow  direction.  Given  the  high  aspect  ratio  and  Peclet 
number  values,  Pe,  it  was  found  reasonable  to  neglect  axial 
dispersion  effects.  Table  1  presents  the  fuel  and  air  channel 
mass  balances.  In  SOFCs,  the  presence  of  CO  and  CO2  in 
the  fuel  does  not  poison  the  anode  reaction,  and  indeed  CO 
can  function  as  a  fuel.  Nevertheless,  here  it  is  assumed  that 
only  hydrogen  is  electrochemically  oxidised  and  that  all  the 
CO  is  converted  through  the  shift  reaction,  considered  to  be 
at  equilibrium  [30-32],  It  is  also  assumed  that  any  methane 
present  in  the  fuel  channel  can  only  be  reformed  to  H2,  CO, 
and  CO2  and,  hence,  not  electrochemically  oxidised.  Thus, 
in  the  fuel  channel,  three  reactions  are  taken  into  account: 
methane  steam  reforming  (i),  water  gas-shift  (ii),  and  hydro¬ 
gen  electrochemical  oxidation  (iii).  In  the  air  channel,  only 
the  reduction  reaction  of  O2  to  O2-  ions  is  considered  (iv). 
Table  2  presents  all  these  reactions.  It  should  be  noted  that 
P(iu)  =  R( iv)  =  A’(v)  where  K(v)  stands  for  the  hydrogen  oxida¬ 
tion  reaction  rate.  Faraday’s  law  relates  the  flux  of  reactants 
and  products  to  the  electric  current  arising  from  an  electro- 
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chemical  reaction.  According  to  this  law,  and  when  only  H2 
oxidation  is  present,  the  local  amount  of  Ho  and  O2  con¬ 
sumed  and  H2O  produced  through  reaction  (v)  in  Table  2  is 
related  to  the  local  electric  current  density,  j,  produced  in  the 
cell  by 

^(v)  =  2F  (1) 

This  relation  relies  on  the  fact  that  solid  electrolytes  such  as 
YSZ  are  pure  ionic  conductors,  and  so  only  oxygen  ions  can 
permeate  through  the  electrolyte  to  take  part  in  the  oxidation 
reaction  at  the  anode/electrolyte  interface  [33]. 

As  for  the  kinetics  of  the  reforming  reaction,  although  con¬ 
ventional  methane  steam  reforming  has  been  widely  studied 
and  there  are  several  publications  dealing  with  its  kinetics 
and  catalytic  aspects  [34—36],  there  is  only  a  small  amount  of 
data  available  for  the  kinetics  of  reforming  in  SOFC  anodes. 
Table  3  presents  some  of  the  rate  equations  published  so  far 
for  Ni  cermet  anodes. 

Lee  et  al.  [12]  studied  methane  conversion  with  differ¬ 
ent  Ni/ZrCL  anode  materials.  For  60  vol.%  Ni  cermets,  they 
found  an  activation  energy  of  98  kJ  mol- 1 .  Achenbach  and 
Riensche  [18]  and  Belyaev  et  al.  [37]  both  showed  that  the 
steam  reforming  reaction  rate  was  first  order  in  methane  and 
zero  order  in  water,  for  a  80  wt.%  ZrC>2  and  20  wt.%  Ni  cer¬ 
met  and  for  a  mixed  Ni-(5  wt.%)  ZrC>2-(2  wt.%)  CeC>2  anode, 
respectively.  Dicks  et  al.  [38]  studied  the  methane  reform¬ 
ing  rate  given  by  a  thin  electrolyte-supported  Ni-YSZ  an¬ 
ode  using  a  tubular  plug  flow  differential  reactor.  Ahmed 
and  Foger  [13]  investigated  Ni-YSZ  and  Ni-YSZ  modi¬ 
fied  (addition  of  a  basic  compound)  anodes.  They  found 
that  the  reaction  order  in  CH4  increased  from  0.85  for  the 


Table  1 

Dynamic  SOFC  model:  mass  balances 
Fuel  channel 


3Cj,f  3C,-,f  1 

ir  =  _Mf^r+  Z.  v‘-^ht 

k6{0),(ii),(v)} 

('  €  {CH4,  h2o,  co,  h2,co2) 

(2) 

cvl=o  =  c°  ie  ich4,  h2o,  co,  h2,  co2) 

(3) 

Air  channel 

dC,-,a  3Qa  D  1  •  rrr  XT  1 

— 7T~  =  Ma  - 1-  Vit(y)R(y)—  l  e  {O2,  N2) 

ot  dx  h  a 

(4) 

Ci, al*=o  =  c?  i  e  {02,  N2) 

(5) 

Table  2 


Reactions  considered  in  the  mass  and  energy  balances 


Steam  reforming  reaction 

n 

+ 

X 

0 

O  CO  +  3H2 

Water  gas- shift  reaction 

(ii)  CO  +  H20  < 

co2  +  h2 

Hydrogen  oxidation  reaction 

(iii)  H2  +  O2-  - 

■>  H20  +  2e" 

Oxygen  reduction  reaction 

(iv)  402  +  2e“ 

-*  O2 

Overall  cell  reaction 

(v)H2  +  ^02  - 

►  h2o 

Table  3 

Rate  equations  for  methane  reforming  on  nickel  cermet  SOFC  anodes 


Kinetic  equation 

Reference 

R0)  =  k0PCH4P^OS  exP 

(6) 

[12] 

(7)%  =  &0PCH4  exp  (-3|r) 

(7) 

[18,37] 

7?(i)  =  ^oPch4  ^1  +  knPi^  +  ks  ^ 

I  exp  (“Sf) 

[38] 

(8) 

(9)S(i)  =  fcoPcH47,H20  exP  (— Mr) 

(9) 

[13] 

Ni-YSZ  anode  to  1.4  for  the  modified  anode,  the  order  in 
steam  changed  from  —0.35  to  —0.8,  and  the  activation  en¬ 
ergy  from  95  to  210kJ  mol-1 .  Other  kinetic  expressions  and 
a  review  in  advances  of  DIR  catalysts  can  be  found  in  refer¬ 
ences  [7,9,1 1,32,39,40].  For  the  fuel  channel  mass  and  energy 
balances,  the  first  order  kinetic  expression  derived  by  Achen¬ 
bach  and  Riensche  [18]  and  presented  in  Table  3,  with  an 
activation  energy  of  Ea  =  82  kJ  mol-1  and  a  pre-exponential 
constant  of  ko  =  4274  mol  s-1  m-2  bar-1,  has  been  adopted. 
Such  a  kinetic  expression  is  considered  typical  of  DIR-SOFC 
performance. 

As  mentioned  previously,  the  water  gas-shift  reaction  (ii) 
is  considered  to  be  at  equilibrium  in  the  fuel  channel;  an 
equilibrium-limited  shift  reaction  rate  expression,  first  or¬ 
der  in  carbon  monoxide  and  with  an  arbitrarily  high  pre¬ 
exponential  factor  Lwgsr.  given  by 

fi(ii)  =  WTO(l-PCOiP"2/PCOP"2°].  00) 
V  ^eq.(ii)  / 

is  used  to  consider  this  effect. 

2.2.  Energy  balances 

In  a  fuel  cell,  current  and  temperature  distributions  are 
strongly  coupled.  Therefore,  knowledge  of  the  prime  heat 
transfer  processes  is  required.  Such  processes  include  heat 
release  arising  from  electrochemical  reactions,  electrical  re¬ 
sistances,  and/or  anode  fuel  chemistry;  convective  heat  trans¬ 
fer  between  cell  components  and  fuel  and  air  gas  streams;  and 
in-plane  heat  conduction  through  cell  components.  Many  au¬ 
thors  have  made  different  assumptions  when  modelling  heat 
transfer  within  a  SOFC.  The  model  in  Table  4  considers  that 
the  thermal  flux  in  the  PEN  and  interconnect  structures  is 
conductive  and  modelled  by  Fourier’s  law  of  heat  conduction, 
and  in  the  gas  channels  is  convective  in  the  gas  flow  direction 
and  from  the  gas  channels  to  the  solid  parts  [1,24].  The  heat 
transfer  coefficients  between  the  gas  channels  and  the  solid 
parts  were  calculated  using  a  constant  Nusselt  number,  Nu 
(a  common  assumption  given  the  laminar  flow  conditions  for 
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which  the  Nusselt  number  is  considered  independent  of  the 
Reynolds  number)  and  the  channel  hydraulic  diameter  [41]. 
Thus,  £f,pEN>  kf  i,  ka.PEN,  and  ka,i  are  given  by: 


^f.PEN  =  £f,i  =  Nu 
^a.PEN  =  £a,I  =  Nil 


<4.f  = 

2Wh( 

(11) 

dh,f  ’ 

W  +  hf 

ka 

^h,a  = 

2  Ml  a 

(12) 

\  a’ 

W  +  ha 

Table  4 

Dynamic  SOFC  model:  energy  balances 


Fuel  channel 


37f  37f  1  1 

PfCp,f-—  =  —UfpfCpf— - b^f,PEN(7pEN  —  Tf)- - b  kfl(T[  —  Tf)  — 

at  ox  h  f  h  f 


-  E 

ks{(i),(ii)l 


(-A  H\Rk 


(13) 


7fh=0  =  if 


(14) 


Air  channel 


97a 


Pacp.a  —  uaPaCp,a  .  P  ^a,PEN(7pEN  Ti)  . 

3f  dx  h a 


+  *a,l(7l  -  7-a)- 
K 


(15) 


7ah=0  =  1? 


(16) 


PEN  structure 

97pen  92rpEN  1 

PPENCp.PEN — - -  =  3.PEN — - Kf,PEN(7pEN  —  Tf) - 

at  ax 1  tpen 


—  ^a,PEN(7pEN  —  7a)- 


1 

rpEN 


+  [(— A77)(V)7?(v)  —  jU] 


a(T] |4  T'pgjj) 


1/ei  +  1/epEN  —  1 


1 

EPEN 

1 

rpEN 


(17) 


where  tpen  —  ranotje  4-  reiectroiyte  T  rcathode 


(IB) 


37pEN 


dx 


=  0, 


37pEN 


dx 


=  0 


(19) 


Interconnect 


3Ti  dzTi  ,  1  1 

PlCp.I—  =  3-1  -r~n - kf,l(T\  —  Tf) - *a,l(7i  “  T-d)  — 

dt  dx~  Ti  Ti 


o(Tf  7^en) 


1/ei +  1/SpEN  -  1 


(20) 


97i 

3a 


=  0, 


37i 

3a 


=  0 


(21) 


where  If  (Xa)  denotes  the  thermal  conductivity  of  the  fuel 
(air)  stream  and  <ih,f  Wh,a)  the  hydraulic  diameter  of  the  fuel 
(air)  channel.  Here,  the  Nusselt  number  is  taken  to  be  3.09 

[42] .  Due  to  the  high  temperatures  involved,  radiation  be¬ 
tween  the  PEN  structure  and  interconnect  is  also  included. 
For  this  purpose,  the  general  expression  for  heat  transfer  by 
radiation  between  two  infinite  parallel  planes  is  considered 

[43] .  It  is  assumed  that  the  hydrogen  oxidation  reaction  en¬ 
thalpy  is  released  at  the  PEN  structure  [29]  and  that,  as  the 
cell  being  modelled  is  placed  in  the  central  region  of  a  stack, 
there  is  no  heat  flux  through  the  external  walls  [22].  Finally, 
it  is  also  assumed  that  the  overall  effect  of  the  variation  of  the 
gas  densities,  heat  capacities  and  velocities  can  be  neglected. 


2.3.  Electrochemical  model 


The  difference  between  the  thermodynamic  potentials  of 
the  electrode  reactions  determines  the  reversible  cell  voltage 
or  open-circuit  potential,  t/°cp  [44,45].  This  open-circuit  po¬ 
tential  is  a  local  quantity,  as  it  depends  on  the  gas  composition 
and  temperature  at  the  electrodes,  and  can  be  determined  by 
the  Nernst  Eq.  (22)  written  for  the  electrochemical  oxida¬ 
tion  of  hydrogen.  stands  for  the  open-circuit  potential  at 
standard  temperature  and  pressure  and  unity  activity,  and  is  a 
function  of  the  Gibbs  free  energy  change  for  the  Ht  oxidation 
reaction. 


t/OCP  =  £/g2 


9tT 

2F 


In 


PH2o.f 

,PH2,fPo22a, 


(22) 


The  theoretical  open-circuit  voltage  is  the  maximum  volt¬ 
age  that  can  be  achieved  by  a  fuel  cell  under  specific  operat¬ 
ing  conditions.  However,  the  voltage  of  an  operating  cell,  U, 
equal  to  the  voltage  difference  between  the  cathode  and  the 
anode,  is  generally  lower  than  this.  As  current  is  drawn  from 
a  fuel  cell,  the  cell  voltage  falls  due  to  internal  resistances 
and  overpotential  losses.  Electrode  overpotential  losses  are 
associated  with  the  electrochemical  reactions  taking  place  at 
the  electrode/electrolyte  interfaces  and  can  be  divided  into 
concentration  and  activation  overpotentials.  These  losses  are 
common  to  all  types  of  fuel  cells  and  cannot  be  eliminated, 
although  temperature,  pressure,  gas  flow  rate  and  composi¬ 
tion,  electrode  and  electrolyte  materials,  and  cell  design,  all 
influence  their  magnitude.  The  voltage  or  potential  of  an  op¬ 
erating  cell  is  thus  generically  given  by  Eq.  (23),  that  can  be 
derived  from  Kirchoff’s  and  Ohm’s  laws  by  considering  a 
cell  as  many  small  elements  in  series,  and  each  one  of  those 
elements  as  a  local  battery  with  its  own  internal  resistance. 

U  =  Uocp  -  (tjoh  m  “I"  ^conc, anode  “1“  ^conc, cathode  H“  ^?act, anode 
"b  *7 act,  cathode)-  (23) 


OOhm  Stands  for  the  ohmic  losses  and  ?7conc, anode?  77conc,  cathode? 
0 act, anode?  and  ^7 act, cathode  for  the  anode  and  cathode  concen¬ 
tration  and  activation  overpotential  losses.  All  these  losses 
are  in  part  responsible  for  the  heat  produced  in  an  operating 
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fuel  cell  and  are  explained  next.  As  both  electrodes  are  nor¬ 
mally  good  conductors,  a  constant  cell  voltage  throughout 
the  cell  is  normally  considered  [46]. 

2.3.1.  Ohmic  losses 

Ohmic  losses  are  caused  by  resistance  to  conduction  of 
ions  (through  the  electrolyte)  and  electrons  (through  the  elec¬ 
trodes  and  current  collectors)  and  by  contact  resistance  be¬ 
tween  cell  components.  This  voltage  drop  is  important  in  all 
types  of  cells  and  is  essentially  linear  and  proportional  to 
current  density.  Because  the  ionic  flow  in  the  electrolyte  and 
the  electronic  flow  in  the  electrodes  obey  Ohm’s  law,  ohmic 
losses  can  be  expressed  by 


Fig.  2.  Schematic  diagram  of  an  anode-supported  solid  oxide  fuel  cell. 


equation  is  written  as: 


POhm  —  j  11  Ohm  -  (24) 

where  j  is  the  current  density  and  Rohm  the  internal  resistance 
of  the  cell,  which  includes  electronic  and  ionic  resistances. 
^Ohm  can  be  obtained  from  measured  cell  resistances  or  be 
estimated  from  the  effective  distance  between  the  compo¬ 
nents  coupled  with  conductivity  data.  Here,  Rohm  is  calcu¬ 
lated  from  the  conductivity  of  the  individual  layers  (assuming 
negligible  contact  resistances,  cross-plane  charge  flow,  and 
series  connection  of  resistances)  by: 

Fanode  ^electrolyte  ^cathode  ,,,, 

Rohm  = - 1 - 1 - ,  (25) 

°anode  ^electrolyte  ^cathode 

where  Tan0de,  teiectrolyte, and  rcathode  stand  for  the  thickness  of 
the  anode,  electrolyte,  and  cathode  layers,  eranode  and  <7Cathode 
for  the  electronic  conductivity  of  the  anode  and  cathode,  and 
o electrolyte  f°r  the  ionic  conductivity  of  the  electrolyte  [47]. 


«  =  <-*§  i" 


PH2O.TPB 


1/2 

,  Ph2,tpb  p  o2,tpb. 


(26) 


The  difference  between  Eqs.  (22)  and  (26)  gives  the  de¬ 
parture  from  the  theoretical  voltage. 


Pconc  — 


mr 

2F 


In 


f  PHoO.TPBPlh.A 
\  PH2  O,  f  Phii ,  TPB  / 


SltF 

4F 


In 


PP2,  a 
Fo2,TPB 


(27) 


where  the  first  term  on  the  right-hand  side  refers  to  the  anodic 
concentration  overpotential,  t?Conc, anode- 


Pconc, anode 


■H 1  ^  (  PH2Q,TPBPH2.f\ 
2  F  V7’H20,fPH2,TPB/ 


(28) 


and  the  second  term  to  the  cathodic  concentration  overpoten¬ 
tial,  Pconc, cathode- 


2.3.2.  Concentration  overpotentials 

Concentration  overpotentials  appear  when  mass  transport 
effects  hinder  the  electrode  reaction,  i.e.  when  the  reactant 
inlet  flux  and  the  product  outlet  flux  from  an  electrode  are 
slower  than  that  corresponding  to  the  discharged  current,  and 
concentration  profiles  develop  across  the  electrodes.  Physical 
processes  that  contribute  to  the  occurrence  of  concentration 
overpotentials  include  gas  species  molecular  transport  in  the 
electrode  pores,  solution  of  reactants  into  the  electrolyte,  dis¬ 
solution  of  products  out  of  the  electrolyte,  and  diffusion  of 
the  reactants/products  through  the  electrolyte  to/from  the  re¬ 
action  sites.  For  an  anode-supported  PEN  structure,  the  con¬ 
centration  overpotential  is  usually  very  small  at  the  cathode 
but  can  be  significant  at  the  anode,  particularly  at  high  current 
densities  and  fuel  utilisation  [48]. 

For  open-circuit  conditions,  the  reactant  and  product  con¬ 
centrations  at  the  electrode/electrolyte  interfaces  or  three- 
phase  boundaries  (TPB)  are  the  same  as  those  in  the  bulk 
channel  flow,  for  which  conditions  the  Nernst  Eq.  (22)  gives 
the  open-circuit  potential.  However,  when  the  current  is  flow¬ 
ing  and  concentration  gradients  develop  (as  illustrated  in 
Fig.  2  for  an  anode-supported  SOFC),  the  species  concen¬ 
trations  at  the  three-phase  boundaries  are  different  from  the 
bulk  concentrations.  In  such  a  case,  the  corresponding  Nernst 


Pconc, cathode 


\  Po2,TPB  / 


(29) 


To  allow  for  the  calculation  of  rjc onc,  a  relation  between 
the  partial  pressures  of  H?,  PbO,  and  CH  at  the  three-phase 
boundaries  and  the  current  density  is  necessary.  Different 
porous-media  gas-phase  transport  models  have  been  de¬ 
veloped  to  predict  concentration  overpotentials  [5,6,48,49]. 
Here,  it  has  been  assumed  that:  the  principal  gaseous  species 
in  the  anode  are  H2  and  bhO  (equimolar  counter-current  one¬ 
dimensional  diffusion  of  H2  and  H2O)  and  in  the  cathode  are 
O2  and  N2  (one-dimensional  self-diffusion);  external  diffu¬ 
sion  is  negligible;  and  the  electrochemical  reaction  occurs 
only  at  the  electrode/electrolyte  interface  [5,6,49].  The  fol¬ 
lowing  equations  were  obtained: 


Ph2,tpb  =  PH2,f  - 


OT7Tanode 
2  F  /-Aff.  anode 


Ph2o,tpb  =  PH20,f  + 


9iTranocje 

2F-Deffanode 


(30) 

(31) 


Po2,tpb  =  P  -  (P  -  Po2,a)  exp 


/  91  TTcathode  A 

y  4 .F Deff ,  cathode  A  / 


(32) 
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where  ran0de  (^cathode)  denotes  the  thickness  of  the  anode 
(cathode),  Z)eff  anotje  represents  the  average  effective  diffusiv- 
ity  coefficient  in  the  anode  (considering  a  binary  gas  mixture 
of  H2  and  H2O),  and  Z)eff, cathode  represents  the  oxygen  effec¬ 
tive  diffusivity  coefficient  in  the  cathode  (binary  gas  mixture 
of  O2  and  N2).  Deff  is  related  to  the  molecular  diffusivity, 
^molecular!  by. 

A:ff,i  =  - - - molecular, i *  *  €  {H2,  H20,  O2,  N2}. 

^tortuosity 

(33) 

This  relation  accounts  for  the  complex  pore  structure  in¬ 
side  the  electrode  material,  where  the  diffusion  path  length 
along  the  pores  is  greater  than  the  measurable  electrode  thick¬ 
ness.  This  is  due  to  the  pores’  tortuous  nature  and  pore  con¬ 
strictions,  and  thus  the  molecular  diffusivity  is  corrected  by  a 
tortuosity  factor,  r tortuosity,  and  the  material  porosity,  ep.  The 
tortuosity  factor  accounts  for  both  the  effects  of  altered  dif¬ 
fusion  path  length  and  changing  cross-sectional  area  in  con¬ 
strictions.  The  molecular  diffusivity  is  the  result  of  fluid-fluid 
intermolecular  collisions  as  considered  by  the  kinetic  theory 
of  gases,  and  is  often  known  as  the  bulk  diffusivity,  D|luik- 
When  the  pore  size  gets  so  small  that  it  is  less  than  the  mean 
path  of  the  fluid,  fluid-wall  collisions  become  the  dominant 
ones  and  the  mode  of  diffusive  transport  is  altered  to  Knud- 
sen  diffusivity,  Z\nudsen,  which  arises  from  the  kinetic  theory 
of  gases  [50].  When  both  these  types  of  diffusion,  bulk  and 
Knudsen,  are  important,  both  contribute  to  the  molecular  dif¬ 
fusivity,  which  is  then  given  by: 

1  1  i 

^molecular,;  T^bulk,/'  f^knudscn./ 

i  e  [H2,  H20,  02,  N2}.  (34) 


2.3.3.  Activation  overpotentials 

Activation  overpotentials  reflect  the  kinetics  of  reactions 
taking  place  on  the  electrode  surface.  At  high  operating  tem¬ 
peratures,  the  electrode  reaction  is  rapid  and,  as  a  result, 
the  activation  overpotential  is  usually  small.  However,  as 
the  operating  temperature  falls,  activation  overpotentials  can 
become  the  most  significant  cause  of  voltage  drop.  Activa¬ 
tion  overpotentials  are  often  represented  by  the  non-linear 
Butler- Volmer  equation,  which  relates  the  current  density 
drawn  to  the  activation  overpotential  and,  for  a  first  order 
charge  transfer  controlled  electrochemical  reaction,  is  given 
by: 


J  —  JO,  electrode 


exp 


anF 
1 WT 


0  act,  electrode 


(1  —  a)nF  ^ 

exp  (  Tj act,  electrode  J 


elec trodee{ anode,  cathode} 


(35) 


where  a  is  the  transfer  coefficient  (usually  taken  to  be  0.5),  n 
the  number  of  electrons  transferred  in  the  single  elementary 


Table  5 

Activation  overpotential  data  for  Eqs.  (35) — (37) 

^cathode  2.35  X  1011  £2_1  m-2  ^cathode  137kJmol_1 

&anode  6.54  X  1011  £2_1  m-2  fianode  140kJmol_l 


rate-limiting  reaction  step  represented  by  the  Butler- Volmer 
equation,  and  70, electrode  the  exchange  current  density,  nor¬ 
mally  expressed  as: 


_  9i  T 

jo, electrode  —  77  ^electrode  CXp 

nr 


^electrode  ^ 


electrode  e  [anode,  cathode}. 


(36) 


However,  in  the  situation  where  charge  transfer  and  mass 
transfer  phenomena  occur  at  comparable  rates,  the  above 
Butler- Volmer  equation  needs  to  be  corrected.  The  extended 
Butler- Volmer  equation  for  the  activation  overpotential  at 
the  anode  (the  cathode  concentration  overpotential  is  subse¬ 
quently  shown  to  be  negligible  and,  thus,  does  not  require 
such  an  extended  approach)  is  then  given  by: 


.1  —  70,  anode 


PH2,TPB 
.  PUi,  f 


exp 


anF 

Inf 


l?act,  anode 


PHeO.TPB 
— - -  exp 

Pn2o,f 


(1  —  a)nF 
9iT~ 


!7act,  anode 


(37) 


For  a  full  knowledge  of  the  activation  overpotentials  in 
a  fuel  cell,  it  is  necessary  to  determine  both  the  exchange 
current  density  and  the  transfer  coefficient  for  both  the 
anode  and  the  cathode  electrode  reactions.  As  previously 
mentioned,  much  research  is  ongoing  to  improve  electrode 
performance.  However,  although  there  are  a  large  number  of 
publications  dealing  with  activation  losses  in  IT-SOFCs,  a 
full  electrode  analysis  that  can  provide  all  the  above  data  in 
a  proven  successful  electrode  is  still  difficult  to  obtain,  as  it 
requires  both  a  detailed  understanding  of  the  electrode  reac¬ 
tion  mechanism  and  its  microstructure.  Here,  for  the  cath¬ 
ode  activation  overpotential,  the  data  provided  by  Horita 
et  al.  [51]  is  used.  This  data  was  obtained  by  studying  the 
cathode  reaction  mechanism  at  the  porous  Lai_  vSrvCo03_d 
(LSC)/Lao.8Sro.2Gao.8Mgo.203_v  (LSGM)  interface  (x  = 
0.2,  0.3,  0.4)  in  the  temperature  range  of  873-1073  K.  This 
represents  one  of  the  few  available  comprehensive  data  sets, 
and  simulations  showed  the  data  to  be  representative  of  more 
conventional  YSZ-based  materials.  The  parameters  n  and  a 
in  Eqs.  (35)— (37)  were  set  equal  to  2  and  0.5,  respectively  (as¬ 
suming  a  two-electron  transfer  rate-determining  step).  After 
fitting  the  available  data  (x  =  0.2),  the  parameters  presented  in 
Table  5  were  obtained.  These  lead  to  a  cathode  area  specific 
resistance  of  0.2  Q  cm2  at  1073  K,  considered  reasonable  for 
IT-SOFC  operation.  For  the  anode,  and  given  that  no  appro¬ 
priate  data  was  found,  it  was  decided  to  use,  as  a  first  estimate, 
the  data  presented  in  Table  5.  The  selection  of  these  data  is 
based  on  the  fact  that  most  authors  indicate  anode  activa¬ 
tion  enthalpies  in  the  range  120-160  kJ  moF* 1  [52,53]  and  in 
the  knowledge  that  an  appropriate  area  specific  resistance  for 
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Table  6 

Model  performance  factors 
Fuel  utilisation  factor 

jLW 

T  Jr  ,  =  _ _ _ 

“  f8^CH4  +  2^2+2Fv°o)Ff0el 


Air  ratio 


v°  F° 

^  _  ^02  ^air 

^•air  —  ="7Trr 


jLW/AF 


3.  Simulation  results  and  discussion 

This  section  is  divided  in  four  sub-sections:  Section  3.1 
introduces  the  model  input  parameters  and  inlet  operating 
conditions;  Section  3.2  presents  the  simulated  curves  for  cell 
voltage  and  power  density  as  a  function  of  current  density; 
Section  3.3  analyses  the  steady-state  performance  of  a  co¬ 
flow  SOFC  system;  and  Section  3.4  briefly  compares  the 
system  steady-state  performance  under  co  and  counter-flow 
operation. 


Power  density 

^sofc  =  jU  (40) 

Fuel  efficiency 


(?CH4^yCH4  +  ALm- 1  +  y°COLHVcoKel 


IT-SOFC  anode  materials  should  not  be  more  than  0. 1  !T2  cm2 
at  1073  K. 


2.3.4.  Performance  factors 

The  resulting  system  of  partial  differential  and  algebraic 
Eqs.  (1)— (6),  (10)— (25),  (28)-(34),  (35)  for  electrode  =  cath¬ 
ode,  and  (36)-(37)  is  solved  using  gPROMS  ModelBuilder 
2.2.5  [54],  with  the  finite  difference  method.  The  initial  con¬ 
ditions  for  the  dynamic  model  (which  is  used  only  on  a  sub¬ 
sequent  publication  to  analyse  the  effect  of  load  changes) 
are  the  initial  steady-state  solution  to  which  a  disturbance 
is  then  imposed.  For  the  solution  of  the  steady-state  model, 
in  addition  to  all  geometry  and  property  data  and  inlet  tem¬ 
peratures  and  compositions,  there  are  three  degrees  of  free¬ 
dom  available  for  the  complete  solution  of  the  system.  One 
option  usually  adopted  in  the  fuel  cell  field  is  to  specify 
the  average  current  density  (or  total  current),  the  fuel  util¬ 
isation  factor,  and  the  air  ratio.  These  are  given  quantities 
from  which  the  molar  fluxes  of  fuel  and  ah'  are  then  de¬ 
termined.  The  fuel  utilisation  factor  is  the  fraction  of  the 
total  inlet  fuel  that  is  used  to  produce  electricity  in  the 
cell.  The  air  ratio  reflects  the  excess  air,  in  relation  to  that 
which  is  stoichiometrically  needed,  that  is  supplied  to  the 
cell  for  cooling.  Output  parameters  usually  adopted  are  the 
power  density  and  the  fuel  efficiency.  The  fuel  efficiency 
represents  the  fraction  of  the  total  chemical  energy  in  the 
inlet  fuel  that  is  converted  into  electrical  energy  (power). 
The  net  chemical  energy  consumed  in  a  fuel  cell  is  con¬ 
verted  either  to  electrical  energy  or  to  heat  generation  due 
to  entropy  change  of  the  electrochemical  reactions,  losses  by 
overpotential  effects,  and  resistance  to  current  flow.  Table  6 
presents  the  mathematical  definition  of  these  performance 
factors. 


3.1.  Input  parameters  and  inlet  operating  conditions 

Table  7  presents  the  model  input  parameters:  geometry 
and  property  data  (properties  of  the  materials  and  cell  di¬ 
mensions),  inlet  temperatures  and  gas  compositions,  average 
current  density  (or  total  current),  fuel  utilisation  factor, 
and  air  ratio.  All  the  physical  property  values  of  the  cell 
materials  are  based  on  mean  values  from  the  literature 
[1,20,22,24,25,42,45].  The  temperature  dependence  of  the 
anode  and  cathode  electronic  conductivity  is  considered  neg¬ 
ligible  in  the  temperature  range  of  concern  [42].  The  PEN 
and  interconnect  emissivity  values  are  considered  typical 
of  ceramic  [42]  and  stainless  steel  materials,  respectively. 
The  cell  dimensions  are  typical  values  for  anode-supported 
intermediate-temperature  SOFCs.  The  inlet  fuel  is  consid¬ 
ered  as  a  gas  mixture  of  CH4,  FEO,  CO,  FE,  and  CO2.  Its 
composition  results  from  a  mixture  with  a  steam  to  carbon 
ratio  equal  to  2  after  10%  pre-reforming,  where  the  shift  re¬ 
action  is  at  equilibrium.  The  average  current  density  value  is 
chosen  based  on  the  curves  presented  in  Section  3.2  and  on 
the  fact  that  fuel  cells  are  normally  designed  to  operate  at  a 
cell  voltage  between  0.6  and  0.7  V.  This  range  of  operation 
is  found  to  be  a  good  compromise  between  cell  efficiency, 
power  density,  low  capital  cost,  stable  operation,  and  avoids 
possible  anode  oxidation  at  low  cell  voltage  [5].  The  air  ra¬ 
tio  value  guarantees  a  temperature  increase  along  the  cell  of 
around  100  K.  The  remaining  model  physical  properties  (gas 
stream  densities,  heat  capacities,  and  thermal  conductivities) 
and  gas  velocities  are  determined  for  the  inlet  conditions,  as¬ 
suming  ideal  gas  behaviour  [42,55],  and  are  taken  as  constant 
values  that  apply  throughout  the  system.  The  thermodynamic 
properties  (reaction  enthalpies  and  Gibbs  free  energy)  are  cal¬ 
culated  as  a  function  of  temperature,  assuming  an  ideal  gas, 
by  integration  (according  to  the  specific  method  rules)  of  the 
heat  capacity  of  each  component  in  the  corresponding  gas 
mixture  [55].  The  effective  diffusivity  coefficients  are  calcu¬ 
lated  assuming  electrodes  with  30%  porosity  (ep  =  0.30),  a 
tortuosity  of  rtortuosity  =  6,  and  an  average  pore  radius  (used 
to  calculate  the  Knudsen  diffusivity)  of  0.5  p,m  [5], 

3.2.  Cell  voltage  and  power  density  as  a  function  of 
current  density 

This  section  presents  the  characteristic  curves  for  cell  volt¬ 
age  and  power  density  as  a  function  of  current  density.  Note 
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Table  7 

Model  input  data  and  operating  conditions 


Physical  properties  of  the  cell  materials 

Anode  electrical  conductivity,  <7anode 

80  x  103  m-1 

Cathode  electrical  conductivity, 

^cathode 

8.4  x  103  Q'1 

PEN  density,  ppen 

5900  kg  m"3 

Interconnect  density,  p\ 

8000  kg  m“3 

PEN  thermal  conductivity,  A.pen 

2  x  10-3  kJ  m-1  s"1  K"1 

Interconnect  thermal  conductivity, 
h 

Interconnect  heat  capacity,  CPii 

25  x  10-3  kJ  m~*  s—  1  K-1 

PEN  heat  capacity,  CPipen 

Electrolyte  ionic  conductivity,  (Xeiectrolyte 

0.5  kj  kg-1  K"1 

33.4  x  103  exp(— 10.3  x  103/T)  ft-1  m-1 

0.5  kJkg~*  K_1 

PEN  emissivity,  Gpen 

0.8 

Interconnect  emissivity,  Gi 

0.1 

Anode  diffusion  coefficient,  Z)eff,  anode 

3.66  x  10-5  m2  s-1 

Cathode  diffusion  coefficient, 

F^eff,  cathode 

1.37  x  10-5  m2  s-1 

Dimensions  of  the  cell  elements 

Cell  length,  L 

0.4  m 

Anode  thickness,  ranode 

500  pm 

Cell  width,  W 

0.1  m 

Cathode  thickness,  rcathode 

50  pm 

Fuel  channel  height,  hf 

1  mm 

Electrolyte  thickness,  reiectroiyte 

20  pm 

Air  channel  height,  ha 

1  mm 

Interconnect  thickness,  i\ 

500  pm 

Operating  conditions 

Fuel  utilisation,  t/fuei 

75% 

Pressure,  P 

1  bar 

Average  current  density,  j 

0.5  Acm~2 

Fuel  inlet  temperature,  Tf° 

1023  K 

Air  ratio,  A.ajr 

8.5 

Air  inlet  temperature,  T® 

1023  K 

Air  feed 

21%  02,  79%  N2 

Fuel  feed 

Steam/carbon  =  2,  10%  pre¬ 
reforming 

that  such  curves  do  not  result  from  the  solution  of  the  full 
SOFC  model  in  Section  2,  but  only  from  the  electrochemical 
model  in  Section  2.3.  The  fuel  composition  in  Figs.  3-5  is 
based  on  a  fully  reformed  steam  and  methane  mixture  (with 
an  initial  steam  to  carbon  ratio  of  2),  to  which  a  degree  of  fuel 
utilisation  (hydrogen  and  carbon  monoxide)  is  then  further 
applied. 

Fig.  3  illustrates  the  predicted  cell  voltage  and  power 
density  as  a  function  of  current  density  at  1073  K  for  an 
undepleted  fully  reformed  fuel  mixture;  the  individual  con¬ 
tributions  of  all  the  various  potential  losses  are  also  shown. 
As  can  be  seen,  the  power  density  reaches  a  peak  value 


of  0.860  W  cm-2  when  the  current  density  is  2.07  A  cm  2 
and  the  cell  voltage  0.415  V.  Under  these  conditions,  the 
cathode  activation  and  ohmic  overpotentials  represent  the 
major  losses  within  the  cell,  followed  by  the  anode  activa¬ 
tion  overpotential.  Although  the  system  under  study  is  an 
anode-supported  one,  anode  concentration  overpotentials  are 
relatively  small,  but  still  much  higher  than  those  at  the  cath¬ 
ode.  A  typical  voltage  versus  current  density  curve  normally 
presents  a  convex  curvature  caused  by  activation  overpoten¬ 
tials  for  low  current  densities  and  a  concave  curvature  asso¬ 
ciated  with  concentration  losses  for  high  current  densities. 
When  there  is  no  fuel  or  oxidant  at  the  three-phase  boundary 
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Fig.  3.  Cell  voltage  and  power  density  as  a  function  of  current  density  at 
1073  K  for  an  undepleted  fully  reformed  fuel  mixture.  The  solid  lines  are 
voltage  curves  (left  axis)  and  the  dashed  line  is  the  power  density  curve  (right 
axis). 


Fig.  4.  Cell  voltage  and  power  density  as  a  function  of  current  density  at 
973,  1023,  and  1073  K  for  an  undepleted  fully  reformed  fuel  mixture.  The 
thinner  lines  are  voltage  curves  (left  axis)  and  the  thicker  lines  are  power 
density  curves  (right  axis). 
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Fig.  5.  Cell  voltage  and  power  density  as  a  function  of  current  density  at 
1073  K  for  a  fully  reformed  fuel  mixture  with  0,  50,  75,  and  85%  fuel  util¬ 
isation.  The  thinner  lines  are  voltage  curves  (left  axis)  and  the  thicker  lines 
are  power  density  curves  (right  axis). 


at  the  electrode/electrolyte  interface,  the  concentration 
overpotentials  become  so  large  that  the  voltage  reaches  zero 
and  the  cell  cannot  operate.  The  current  density  at  which  this 
phenomenon  occurs  is  often  referred  to  as  the  limiting  current 
density  (anode  or  cathode  limiting  current  density).  However, 
whether  or  not  a  concave  curvature  is  observed  at  high  current 
densities  depends  on  the  particular  cell  operating  conditions, 
design,  and  materials,  which  determine  the  relative  magni¬ 
tude  of  the  ohmic  and  activation  losses.  For  the  cell  and  op¬ 
erating  conditions  illustrated  in  Fig.  3,  no  concave  curvature 
is  observed.  This  is  mainly  due  to  high  ohmic  and  activation 
losses  that  cause  the  cell  voltage  to  drop  to  zero  when  the  cur¬ 
rent  density  reaches  4.45  A  cm-2  (while  the  anode  limiting 
current  density  of  10.1  A  cm-2  is  never  reached).  As  men¬ 
tioned  previously,  fuel  cells  are  normally  designed  to  operate 
at  a  cell  voltage  between  0.6  and  0.7  V.  Thus,  if  the  cell  under 
study  were  to  be  operated  at  a  constant  temperature  and  gas 
composition,  the  best  operating  point,  based  on  Fig.  3,  would 
be  between  current  density  values  of  0.823  and  1 .23  A  cm-2, 
which  lead  to  power  densities  between  0.583  and 
0.739  W  cm-2. 

All  the  above  is  valid  for  operation  at  1073  K.  However, 
and  as  pointed  out  before,  one  of  the  main  incentives  in  us¬ 
ing  anode-supported  SOFC  designs  is  to  allow  a  reduction 
in  the  operating  temperature  and,  thus,  the  use  of  both  metal 
and  ceramic  materials  as  compared  to  all-ceramic  SOFCs. 
Thus,  it  is  important  to  evaluate  the  cell  performance  at  lower 
temperatures.  Fig.  4  illustrates  the  predicted  cell  voltage  and 
power  density  as  a  function  of  current  density  at  973,  1023, 
and  1073  K,  for  the  same  fuel  composition  as  in  Fig.  3.  As 
expected,  the  performance  of  a  cell  is  significantly  hindered 
when  the  operating  temperature  is  decreased:  the  maximum 
power  density  reduces  from  0.860  W  cm-2  at  1073  K,  to 
0.539  W  cm-2  at  1023  K  (at  a  current  density  of  1.36  A  cm-2 
and  a  voltage  of  0.397  V)  and  to  0.314  W  cm-2  at  973  K  (at 


a  current  density  of  0.823  A  cm-2  and  a  voltage  of  0.381  V); 
and  the  point  at  which  the  cell  voltage  drops  to  zero  reduces 
from  4.45  A  cm-2  at  1073  K  to  2.93  A  cm-2  at  1023  K  and 
to  1.79  A  cm-2  at  973  K.  The  best  operating  point,  applying 
the  same  criteria  as  above,  now  falls  between  current  density 
values  of  0.518  and  0.734Acm-2  for  1023 K  (which  leads 
to  power  densities  between  0.360  and  0.443  W  cm-2)  and 
between  0.269  and  0.428  A  cm-2  for  973  K  (which  leads  to 
power  densities  between  0.190  and  0.253  W  cm-2). 

Fuel  composition  plays  a  significant  role  in  fuel  cell  per¬ 
formance.  In  operation  and  as  the  fuel  flows  through  a  SOFC 
fuel  channel,  hydrogen  and  carbon  monoxide  are  consumed 
by  the  electrochemical  reactions  and  the  fuel  stream  be¬ 
comes  diluted,  being  less  rich  close  to  the  exit  of  the  chan¬ 
nel.  Therefore,  in  this  region,  concentration  overpotentials 
become  more  significant.  Fig.  5  illustrates  the  predicted  cell 
voltage  and  power  density  as  a  function  of  current  density 
for  a  fully  reformed  fuel  mixture  with  0,  50,  75,  and  85% 
fuel  utilisation.  For  0  and  50%  fuel  utilisation,  concentra¬ 
tion  overpotentials  are  not  significant.  As  seen  in  Fig.  5, 
for  these  two  cases,  the  voltage  versus  current  density  curve 
does  not  present  a  concave  curvature  for  high  current  den¬ 
sities.  However,  for  75  and  85%  fuel  utilisation,  the  con¬ 
cavity  becomes  visible  and  voltage  and  power  density  drop 
more  rapidly  to  zero,  due  to  concentration  overpotentials 
in  the  anode  (cathode  concentration  overpotentials  remain 
negligible).  The  anode  limiting  current  density  decreases 
significantly  from  5.64  A  cm-2  at  50%  fuel  utilisation,  to 
2.98  A  cm-2  at  75%  and  1.83  A  cm-2  at  85%.  However,  as 
the  activation  overpotentials  also  take  into  account  the  effect 
of  mass  transport  through  Eq.  (37),  the  point  at  which  the 
cell  voltage  drops  to  zero  occurs  before  (at  3.59  A  cm-2  for 
50%  fuel  utilisation,  2.55  A  cm-2  for  75%  and  1.71  A  cm-2 
for  85%). 

In  terms  of  reported  values  within  the  available  literature, 
main  results  are:  peak  power  densities  of  1.08  W  cm-2  at 
973  K  and  1.58  W  cm-2  at  1073  K  [49]  or  0.49  W  cm-2  at 
1073  K  [5]  for  anode-supported  SOFCs  operating  on  humid¬ 
ified  hydrogen,  and  0.40  W  cm-2  at  1023  K  for  an  anode- 
supported  SOFC  operating  on  methane  [48].  Important  to 
note  is  the  fact  that,  although  very  high  power  densities  have 
been  reported  for  single  cell  test  by  SOFC  developers,  no 
data  is  yet  available  to  allow  for  the  validation  of  developed 
models.  In  addition,  it  should  also  be  noted  that  power  den¬ 
sity  values  for  stack  operation  are,  as  confirmed  in  Section 

3.3,  well  below  single  cell  values. 

3.3.  Anode-supported  IT  DIR-SOFC  operation 

The  model  described  in  Section  2  is  able  to  predict  the  var¬ 
ious  temperatures  along  the  cell  length  (fuel  and  air  channels, 
PEN  structure,  and  interconnect),  the  gas  composition  in  the 
fuel  and  air  channels,  all  the  electrochemical-related  variables 
(open-circuit  voltage,  activation,  ohmic,  and  concentration 
overpotential  losses,  terminal  potential,  and  current  density) 
as  well  as  the  cell  efficiency  and  power  output.  Figs.  6-8 
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Fig.  6.  Fuel  channel  component  mole  fractions  along  the  cell  length  for  the  Fig-  8.  Predicted  voltage,  current  density,  and  contribution  of  all  the  various 

conditions  in  Table  7.  potential  losses  along  the  cell  length  for  the  conditions  in  Table  7. 


present  some  of  these  variables  for  the  operating  conditions 
given  in  Table  7.  Figs.  6-8  are  for  a  co-flow  operation  case, 
while  Figs.  12-14,  in  Section  3.4,  are  for  counter-flow  oper¬ 
ation. 

Fig.  6  presents  the  mole  fraction  profiles  along  the  cell 
length  of  all  the  components  in  the  fuel  channel  stream.  These 
illustrate  the  impact  of  the  simultaneous  occurrence  of  the  di¬ 
rect  internal  reforming  reaction,  the  water  gas-shift  reaction, 
and  the  oxidation  of  hydrogen  at  the  anode/electrolyte  inter¬ 
face.  At  the  entrance  of  the  fuel  cell,  and  due  to  the  high 
methane  content,  the  reforming  reaction  is  much  faster  and 
the  methane  starts  being  consumed  rapidly,  producing  hydro¬ 
gen  and  carbon  monoxide.  Once  most  of  the  methane  is  con¬ 
sumed,  the  hydrogen  oxidation  reaction  becomes  the  faster 
reaction  and  the  consumption  of  hydrogen  and  production  of 
steam  can  be  clearly  seen.  The  current  density  profile  along 
the  cell  length,  shown  in  Fig.  8,  corroborates  the  increasing 


Fig.  7.  Fuel  and  air  channels,  PEN  structure,  and  interconnect  temperature 
along  the  cell  length  for  the  conditions  in  Table  7. 


hydrogen  oxidation  reaction  along  the  cell.  At  the  exit  of  the 
fuel  channel,  all  the  methane  has  been  fully  consumed  and 
the  stream  content  is  65%  in  FbO,  4%  in  CO,  15%  in  FF,  and 
16%  in  C02. 

Fig.  7  presents  the  fuel  and  air  channels,  PEN  structure, 
and  interconnect  temperature  profiles  along  the  cell  length.  It 
can  be  seen  that  the  cell  temperature  increases  along  the  fuel 
and  air  flow  directions  by  heat  accumulation,  with  the  max¬ 
imum  temperature  occurring  at  the  outlet.  This  temperature 
distribution  is  typical  of  a  co-flow  arrangement.  Although,  the 
steam  reforming  reaction  is  known  to  be  strongly  endother¬ 
mic,  it  can  be  seen  that  there  is  enough  heat  production  in 
the  cell  to  enable  the  reaction.  The  total  heat  consumed  by 
the  DIR  reaction  represents  45%  of  the  total  heat  produced 
in  the  cell.  The  simulations  do  illustrate  a  slight  local  cooling 
effect  associated  with  the  reforming  of  methane,  where  the 
fuel  channel  temperature  decreases  by  54  K  in  relation  to  its 
inlet  value  (the  equivalent  temperature  decrease  in  the  PEN 
structure  is  only  18  K).  As  mentioned  before,  direct  or  indi¬ 
rect  internal  reforming  in  typical  high-temperature  SOFCs, 
where  the  inlet  temperature  is  normally  1173K,  is  charac¬ 
terised  by  severe  local  cooling  effects  at  the  cell  entrance. 
The  profiles  in  Fig.  7  then  demonstrate  that  intermediate- 
temperature  SOFCs  can  represent  a  great  improvement  in  this 
area.  Also,  note  that  the  extent  of  pre-reforming  assumed  for 
these  simulations  is  of  only  10%,  and  thus  the  inlet  fuel  is  very 
rich  in  methane.  The  maximum  and  minimum  temperatures 
in  the  PEN  structure  are  1135  and  981  K,  respectively,  lead¬ 
ing  to  a  maximum  temperature  difference  of  154  K  in  the 
cell.  Temperature  gradients  in  a  solid  oxide  fuel  cell  stack 
(along  the  length  and  across  the  various  layers)  need  to  be 
carefully  monitored  to  avoid  any  thermally  induced  fractures 
on  its  ceramic  components,  with  a  consequent  cell  failure. 

Finally,  Fig.  8  shows  the  predicted  voltage  and  current 
density  along  the  cell  length,  as  well  as  the  individual  con¬ 
tribution  of  all  the  various  potential  losses.  As  discussed  be¬ 
fore,  for  a  given  cell,  the  potential  losses  and  current  density 
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(b)  Dimensionless  axial  position 


Fig.  9.  Fuel  and  air  channels,  PEN  structure,  and  interconnect  temperature  along  the  cell  length  for  the  conditions  in  Table  7,  except  that  the  fuel  and  air  inlet 
temperatures  are:  (a)  973  K  and  (b)  1073  K. 


distribution  along  the  cell  are  related  to  the  composition  of  the 
gas  streams  and  to  the  various  temperatures  in  the  cell.  The 
open-circuit  potential  is  governed  by  the  Nernst  equation,  and 
thus  follows  both  the  cell  temperature  and  the  partial  pressure 
of  hydrogen,  as  this  is  produced  by  DIR  and  then  oxidised 
and  replaced  by  steam.  As  the  temperature  increases  along 
the  cell  length,  the  overpotentials  decrease,  allowing  for  an 
increase  of  the  current  density.  The  concentration  overpoten¬ 
tials,  even  though  very  small  in  the  case  under  study,  increase 
as  the  fuel  is  consumed  along  the  cell.  In  percentage  terms, 
the  cathode  activation  overpotentials  represent  45%  of  the 
total  losses  within  the  cell,  followed  by  the  anode  activation 
overpotentials  with  31%,  and  the  PEN  structure  ohmic  losses 
with  22%.  The  anode  concentration  overpotentials  contribute 
to  2%  of  the  losses  while  the  losses  by  cathode  concentra¬ 
tion  overpotentials  were  found  negligible  (expected,  given 
the  reduced  cathode  thickness  and  oxygen  utilisation).  This 
is  the  expected  trend  for  an  IT-SOFC,  where  the  electrolyte 
is  thin  and  ohmic  losses  no  longer  dominate.  For  the  op¬ 


erating  conditions  in  Table  7,  the  SOFC  under  study  has  a 
satisfactory  terminal  voltage  of  0.663  V,  a  power  density  of 
0.332  W  cm-2,  and  a  fuel  efficiency  of  46.8%. 

Figs.  9-11  illustrate  the  effect  of  changing  some  of  the 
operating  conditions.  Fig.  9  presents  the  fuel  and  air  chan¬ 
nels,  PEN  structure,  and  interconnect  temperature  profiles 
along  the  cell  length  for  the  cases  where  the  fuel  and  air  inlet 
temperatures  are  973  and  1073  K  (all  the  remaining  condi¬ 
tions  in  Table  7  are  kept  constant).  From  Figs.  7  and  9,  it 
can  be  seen  that  as  the  inlet  temperature  increases,  the  tem¬ 
perature  gradients  at  the  entrance  of  the  fuel  channel  also 
increase.  This  is  due  to  an  increase  in  the  steam  reforming 
reaction  rate  with  higher  temperatures.  However,  the  overall 
temperature  increase  across  the  cell  length  decreases.  This  is 
because,  at  a  given  current  density,  higher  temperatures  lead 
to  a  more  electrically  efficient  operation  (the  fuel  cell  effi¬ 
ciency  is  40.2%  for  fuel  and  air  inlet  temperatures  of  973  K 
and  52.0%  for  1073  K),  and  thus  to  a  lower  heat  generation 
in  the  cell.  This  implies  that  if  the  studied  SOFC  needs  to 
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Fig.  10.  Predicted  voltage,  current  density,  and  contribution  of  all  the  various  potential  losses  along  the  cell  length  for  the  conditions  in  Table  7,  except  that  the 
fuel  utilisation  is:  (a)  85%  and  (b)  95%. 
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(a)  Dimensionless  axial  position 
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Fig.  11.  Fuel  and  air  channels,  PEN  structure,  and  interconnect  temperature  and  predicted  voltage,  current  density,  and  individual  contribution  of  all  the  various 
potential  losses  along  the  cell  length  for  the  conditions  in  Table  7,  except  that  the  average  current  density  is  0.75  A  cm  2. 


be  operated  at  lower  temperatures  while  keeping  the  overall 
temperature  gradient  across  the  cell  of  100K,  the  air  ratio 
needs  to  be  increased. 

Fig.  10  presents  the  predicted  voltage,  current  density,  and 
individual  contributions  of  all  the  various  potential  losses 
along  the  cell  length  for  the  cases  where  the  fuel  utilisation  is 
changed  to  85  and  95%.  The  purpose  of  showing  Fig.  10  is  to 
illustrate  the  effect  of  concentration  overpotentials  on  the  cell 
performance.  As  can  be  seen,  as  the  fuel  is  consumed  along 
the  cell,  concentration  overpotentials  become  more  signifi¬ 
cant  and  the  current  density  that  can  be  locally  drawn  drops 
considerably  to  maintain  a  constant  voltage.  This  effect  is 
also  accentuated  because  the  open-circuit  voltage  decreases 
significantly  with  the  consumption  of  hydrogen.  As  expected, 
the  fuel  efficiency  increases  with  an  increase  in  fuel  utilisa¬ 
tion  (53.6%  for  a  fuel  utilisation  of  85%  and  59.8  for  95%). 
However,  and  since  the  amount  of  inlet  fuel,  and  thus  the 
amount  of  heat  consumed  to  reform  all  the  methane  is  much 
lower,  the  overall  temperature  profiles  are  similar  to  those 
presented  in  Fig.  7. 

Fig.  1 1  presents  the  fuel  and  air  channels,  PEN  structure, 
and  interconnect  temperature  and  predicted  voltage,  current 
density,  and  individual  contribution  of  all  the  various  poten¬ 
tial  losses  along  the  cell  length  for  the  case  where  the  cell  is 
operated  at  an  average  current  density  of  0.75  A  cm-2.  It  can 
be  seen  that,  with  an  increase  in  current  density,  the  cell  tem¬ 
perature  increases.  This  implies  that  to  maintain  a  constant 
total  temperature  increase  along  the  cell,  an  increase  in  the  air 
ratio  is  required.  From  Fig.  1  lb,  it  can  be  seen  that  an  increase 
in  current  density  leads  to  an  increase  in  all  the  overpoten¬ 
tial  losses  in  a  cell  and  thus  to  a  lower  operating  voltage.  At 
0.75  A  cm-2,  the  cell  voltage  decreases  to  0.606  V,  the  power 
density  is  0.455  W  cm-2,  and  the  fuel  efficiency  42.7%.  In 
summary,  temperatures,  fuel  composition,  gas  flowrates,  and 
current  density,  greatly  affect  the  cell  performance  and  source 
of  operation  losses.  The  best  combination  of  operating  con¬ 
ditions  will  always  be  a  compromise  of  several  factors,  not 


only  related  to  the  cell  stack  but  to  all  the  balance  of  plant 
involving  it. 

3.4.  Anode-supported  IT  DIR-SOFC  operation:  co-flow 
versus  counter-flow  operation 

This  last  section  aims  to  present  a  brief  comparison  of 
the  system  behaviour  under  co  and  counter-flow  operation. 
The  air  channel  mass  and  energy  balances  and  corresponding 
boundary  conditions  (Eqs.  (4)  and  (5)  in  Table  1  andEqs.  (15) 
and  (16)  in  Table  4)  have  been  modified  accordingly  and  the 
same  operating  conditions  as  the  ones  given  in  Table  7  have 
been  considered. 

Figs.  12-14  present  the  fuel  and  air  channels,  PEN  struc¬ 
ture,  and  interconnect  temperature  profiles,  the  mole  fraction 
profiles  of  all  the  components  in  the  fuel  channel  gas  stream, 
and  the  voltage,  current  density,  and  individual  contribution 
of  all  the  various  potential  losses,  along  the  cell  length  for 


Fig.  12.  Fuel  and  air  channels,  PEN  structure,  and  interconnect  temperature 
along  the  cell  length  for  the  conditions  in  Table  7  for  counter-flow  operation. 
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Fig.  13.  Fuel  channel  component  mole  fractions  along  the  cell  length  for  the 
conditions  in  Table  7  for  counter-flow  operation. 


counter-flow  operation.  Figs.  6-8  above  present  the  equiva¬ 
lent  co-flow  operation  figures. 

Comparing  Figs.  7  and  12  it  can  be  seen  that,  while  in  a 
co-flow  SOFC  the  temperature  increases  along  the  air  and 
fuel  flow  direction  by  heat  accumulation,  with  the  maximum 
temperature  at  the  exit,  in  a  counter-flow  cell  the  temper¬ 
ature  profile  has  a  peak  near  the  air  channel  exit  and  the 
operation  is  characterised  by  higher  average  cell  tempera¬ 
tures  and  steeper  local  temperature  gradients.  Comparing 
Figs.  6  and  13  it  is  observed  that,  under  counter-flow  op¬ 
eration,  and  since  the  temperature  close  to  the  fuel  channel 
entrance  is  higher  (Fig.  12),  all  the  methane  is  fully  reformed 
in  the  first  20%  of  the  cell  length,  while  in  the  co-flow  case, 
methane  is  gradually  reformed  along  the  fuel  channel  length. 
As  discussed  above,  temperature  gradients  in  a  SOFC  stack 
need  to  be  carefully  monitored  to  avoid  any  thermally  induced 
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Fig.  14.  Predicted  voltage,  current  density,  and  contribution  of  all  the  various 
voltage  losses  along  the  cell  length  for  the  conditions  in  Table  7  for  counter¬ 
flow  operation. 


fractures  on  its  ceramic  components.  This  seems  to  represent 
the  major  drawback  of  SOFC  operation  under  counter-flow 
of  the  fuel  and  air  gas  streams:  in  the  co-flow  operation  case 
the  simulations  illustrate  only  a  slight  local  cooling  effect  as¬ 
sociated  with  DIR  (Fig.  7);  but  in  the  counter-flow  case  the 
temperature  has  a  steep  135  K  increase  in  the  first  8  cm  of 
the  fuel  channel.  Fig.  14  illustrates  once  again  how  the  var¬ 
ious  voltage  losses  and  the  current  density  distribution  are 
related  to  gas  composition  and  cell  temperature:  the  open- 
circuit  voltage  follows  both  the  cell  temperature  and  partial 
pressures  of  hydrogen  and  steam;  and  the  concentration  over¬ 
potentials,  even  though  negligible  for  the  simulated  operat¬ 
ing  conditions,  are  more  pronounced  for  higher  local  current 
densities. 

Finally,  for  the  operating  conditions  in  Table  7,  this 
counter-flow  SOFC  system  has  a  terminal  voltage  of  0.771  V, 
a  power  density  of  0.386  W  cm-2,  and  a  fuel  efficiency  of 
54.4%.  The  considerably  better  performance  observed  for 
the  SOFC  system  operating  under  counter-flow  is  an  out¬ 
come  of  the  higher  average  cell  temperatures  and  consequent 
lower  cell  activation  overpotentials  and  ohmic  losses.  How¬ 
ever,  because  of  the  steep  temperature  gradients  in  Fig.  12, 
not  desirable  for  a  ceramic  system,  and  the  consequent  un¬ 
even  current  density  distribution  (Fig.  14),  counter-flow  is 
not  considered  the  best  mode  of  operation.  Achenbach  [29] 
reached  similar  conclusions  for  a  high-temperature  SOFC 
system. 


4.  Conclusions 

A  dynamic  anode-supported  intermediate-temperature  di¬ 
rect  internal  reforming  planar  one-dimensional  SOFC  stack 
model  has  been  presented.  The  model  developed  consists 
of  mass  and  energy  balances,  and  an  electrochemical  model 
that  relates  the  fuel  and  air  gas  composition  and  temperature 
to  voltage,  current  density,  and  other  relevant  fuel  cell  vari¬ 
ables.  The  electrochemical  performance  of  the  cell  has  been 
analysed  for  several  temperatures  and  fuel  utilisations.  Max¬ 
imum  power  densities  of  0.86  W  cm-2  for  a  current  density 
of  2. 1  A  cm-2  and  a  cell  voltage  0.42  V  have  been  reported 
for  an  undepleted  fully  reformed  fuel  mixture  at  1073  K.  The 
steady-state  performance  of  the  cell  as  well  as  the  impact  of 
changes  in  inlet  temperatures,  fuel  utilisation,  average  cur¬ 
rent  density,  and  flow  configuration  have  been  illustrated. 
The  SOFC  system  studied  operates  at  an  output  voltage  of 
0.66  V,  a  power  density  of  0.33  W  cm-2,  and  a  fuel  efficiency 
of  47%,  when  using  a  10%  pre-reformed  methane  fuel  mix¬ 
ture  with  75%  fuel  utilisation,  inlet  fuel  and  air  temperatures 
of  1023  K,  average  current  density  of  0.5  A  cm-2,  and  an  air 
ratio  of  8.5.  It  has  been  shown  that,  for  similar  operating  con¬ 
ditions,  SOFC  operation  under  counter-flow  of  the  fuel  and  air 
gas  streams  leads  to  steep  temperature  gradients  with  a  con¬ 
sequent  uneven  current  density  distribution,  and  hence,  is  not 
considered  the  best  mode  of  operation.  The  model  presented 
is  used,  in  a  subsequent  publication,  to  analyse  the  dynamic 
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response  of  the  described  anode-supported  IT  DIR-SOFC  to 
possible  load  changes  [19]. 
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